Slippage of water past superhydrophobic carbon nanotube forests in microchannels 
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We present in this letter an experimental characterization of Uquid flow slippage over superhy- 
drophobic surfaces made of carbon nanotube forests, incorporated in microchannels. We make use 
of a /i-PIV (Particule Image Velocimetry) technique to achieve the submicrometric resolution on 
the flow profile necessary for accurate measurement of the surface hydrodynamic properties. We 
demonstrate boundary slippage on the Cassie superhydrophobic state, associated with slip lengths 
of a few microns, while a vanishing slip length is found in the Wenzel state, when the liquid im- 
pregnates the surface. Varying the lateral roughness scale L of our carbon nanotube forest-based 
superhydrophobic surfaces, we demonstrate that the slip length varies linearly with L in line with 
theoretical predictions for slippage on patterned surfaces. 
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The nature of the hydrodynamic boundary condition 
of a liquid past a solid surface has been actively revisited 
over the recent years, strongly motivated by the emer- 
gence of microfluidics and nanofluidics 0, Q • Downsiz- 
ing leads to large surface-to- volume ratio, and the flow 
properties at the solid boundaries become a key factor in 
the understanding of the motion of liquids in ever smaller 
channels In this context the violation of the no-slip 
boundary condition (BC) and the existence of slippage of 
the fluid at the solid surface is a promising way to bypass 
the huge increase in hydrodynamic resistance that comes 
with miniaturization. 

Slippage, usually quantified by a slip length 6, has been 
demonstrated on solvophobic surfaces^with a slip length 
of a few tens of nanometers |2fc S S • Such values 
are obviously insufficient to modify the flow in channels 
with micrometer sizes and other specific solutions have 
to be proposed to benefit from slippage effects in mi- 
crofluidics. Super- hydrophobic surfaces are promising in 
this context, since it has been recently predicted and 
experimentally reported 8] that the associated "Fakir" 
effect (the so-called Cassie state) considerably amplifies 
slippage. This situation is achieved with highly rough 
hydrophobic surfaces: instead of entering this bidimen- 
sionnal hydrophobic porous medium, the liquid remains 
at the top of the roughness thus trapping some air in the 
interstices, and therefore leading to a very small liquid- 
solid contact area. 

In the context of microfluidics, the challenge now con- 
sists in developping versatile methods to design such sur- 
faces in microchannels with optimized flow properties. 
These have to deal with two conflicting constraints on 
the engineered surfaces: low friction -i.e. large slippage- 
and robustness of the Fakir effect against pressure in- 
duced impregnation. On the one hand, large slippage 
is achieved by minimizing the liquid-solid area together 



with maximizing the underlying lateral length scale L of 
the roughness. Indeed, theoretical predictions show that 
the resulting effective slip len gth 6pff is mainly fixed by 
the roughness scale L j^l, lid UjJ, ^eff — with 
a a numerical factor (smaller than unity for typical sur- 
faces). On the other hand, the stability of the Fakir 
(Cassie) state with respect to pressure variations fixes 
some restrictive conditions on the underlying roughness 
scale L 0: above an excess pressure of the order of 
APint ^ 2jly/L (with the liquid- vapor surface ten- 
sion), penetration of the liquid into the roughness occurs 
-Wenzel state-, and the benefits from the low friction 
Fakir state are lost. The potential applicability of litho- 
graphically designed patterned surfaces, with a surface 
lateral scale L in the ten microns range is therefore lim- 
ited by their low resistance to pressure j^. While one 
would expect the increased robustness on nanoengineered 
surfaces to be achieved at the expense of their slippage 
properties, recent rheology measurements on nanopost 
superhydrophobic surfaces have reported slip lengths of 
tens of microns [iJ in clear contrast with theoretical ex- 
pectations 1^, lid IllL and with the experimental re- 
sults of at larger (micrometric) scales. This surprising 
result would suggest that the underlying roughness does 
not limit the slip length for nanopatterned surfaces, as 
first expected. However the sensitivity of the rheologi- 
cal method used in 14] appears questionnable and 
it therefore points out the need for accurate experimen- 
tal data to achieve a proper understanding of the role 
of the underlying roughness scale, from nano- to micro- 
patterned surfaces. 

The purpose of the present work is precisely to pro- 
vide accurate data allowing to understand the role of the 
roughness scale on the hydrodynamic properties of su- 
perhydrophobic surfaces. Taking benefit of recent ad- 
vances in nanomaterials science we incorporated ver- 
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FIG. 1: Experimental setup. A pressure driven flow is con- 
ducted in a glass/NTC micro- channel. The focal plane Z is 
controlled with a piezo, a large numerical aperture objective 
allows for a narrow depth of field. Velocity is measured by 
particle image velocimetry over the whole channel depth. 

satile composite surfaces in microchanels. We make 
use of recently developed surfaces covered with verti- 
cally ali gne d and densely spaced carbon nanotubes (CNT 
forests) [lall3- Together with the existence of large con- 
tact angles (larger than 165°), the submicrometric lateral 
length scale achievable on these surfaces was shown to 
provide a robust Fakir state with respect to pressure ^3 ^ 
with an intrusion pressure Pint ^ latm. More interest- 
ing in the present framework, is the ability provided by 
the method to achieve various lateral length scales L in 
the submicron to microns range. Flows in the vicinity of 
such surfaces were then investigated using particle image 
velocimetry at the micron scale (/i-PIV) to probe the ve- 
locity profile in the microchannels with submicrometric 
resolution. Our measurements confirm that a low fric- 
tion is achieved providing the interface lies in the Fakir 
state, and demonstrate that these low friction properties 
are indeed controlled by the lateral roughness scale, in 
agreement with theoretical considerations 0, . 

The experimental setup is shown in fig. The CNT 
forests are obtained using Plasma Enhanced Chemical 
Vapor Deposition (PECVD), using the protocol previ- 
ously described in 0, ^3 • The CNT forests are grown 
on a silicon substrate covered with a Si02 layer and a thin 
film of Ni catalyst. This develops Ni islands through the 
sintering of the Ni film (after introduction in a furnace 
heated to 750° C under 20 torr of hydrogen). Ammo- 
nia and then acethylene (C2H2) are introduced into the 
furnace to initiate the CNT growth. A dc discharge be- 
tween the cathode and anode samples allows for directed 
growth. The depositions are carried out between 5 and 
60 minutes in a stable discharge. Typical CNT forests 
grown through this process are characterised by a typical 
nanotube diameter between 50 and 100 nm, a nanotube 
inter-distance 100-250 nm, and micrometric lengths, as 
shown on Fig. Efa). The nanotube lengths and inter- 
distance can be tuned by varying the deposition time 
and the thickness of the nickel layer. 




FIG. 2: SEM pictures of superhydrophobic CNT forests, after 
functionalization with thiols: (a) in the gas phase; (b-d) in 
ethanol, for increasing initial nanotube length, resulting in 
characteristic lateral roughness length scale L: 1.7, 3.5 and 
6 /xm (pictures b, c and d). 

To get a superhydrophobic behavior, the CNT forests 
are then functionalized using thiols in gas or liquid phase, 
after covering them with an anchoring thin gold layer (5 
nm) by sputtering as described in 17]. While the gas 
phase functionalization leaves the CNT forest structure 
unchanged 17], the functionalization in a liquid phase 
results in bundles of carbon nanotubes (see figl^^b-d)), 
as a result from capillary adhesion during the eva por ative 
drying of the functionalization solvent, ethanol [iTl Il8j. 
Combined with the control of the nanotube density and 
length, this bundling process is used in the present study 
as a further way to tune the roughness lateral scale L. 

These CNT surfaces are embedded in a microfiuidic 
setup (see fig. ^ featuring a microchannel of height 40- 
50/im, width 2mm and length 2cm. This microchannel 
is obtained by sandwiching a dry film resist that has a 
lithographic trench between a CNT covered silicon sub- 
strate and a microscope glass coverslip, following a re- 
cently introduced method A pressure driven flow 
is set up by controlling inlet and outlet pressures in the 
range of -50 to 100 mbar with respect to the atmospheric 
pressure, while the pressure drop along the channel re- 
mains below 2 mbar. The working solution is deionized 
water seeded with fluorescent particles (500 nm in diam- 
eter), 10~^ in volumetric concentration. The observation 
is made with a LEICA epifluorescent microscope, using 
a 100 X oil immersion large numerical aperture objective 
(NA=1.3) achieving a 700nm depth of field. An optical 
depth of 500nm is obtained for the measurement volume 
by thresholding on the intensity and then selecting the 
particles in focus. The velocity profile is measured using 
a standard /i-PIV technique, similar to that described in 
[2^, where further details can be found. Briefly, a scan 
in the vertical z direction is provided using a piezotrans- 
ducer with 0.2 jam step increments, allowing to measure 
the velocity profile v{z) in the whole microchannel. A 
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time- delayed intensity cross correlation is performed in a 
20/imx60/im window, and the computed velocity is av- 
eraged over 20 pairs of images for each z-position of the 
piezo. We thus obtain eventually an averaged velocity in 
an imaged zone of 500nm x 20/im x 60/im. Typical veloc- 
ity profiles (see fig. ^ present the characteristic parabolic 
shape associated with the pressure-driven Poiseuille flow. 

The position of the CNT surface has been carefully 
determined. As in [2^, we take advantage of the pres- 
ence of a few particles which unavoidably adsorb on sur- 
faces [we only considered windows of observation where 
at least 3, and up to ten, of such beads are included]. 
The position of the surface is then determined by fitting 
the fluorescence intensity emitted by these beads with a 
lorentzian and by substracting the bead radius from the 
location of the maximum of the lorentzian. We have con- 
sidered flows on surfaces in the Cassie state (with trapped 
air), but also on surfaces in the Wenzel state, where the 
CNT surfaces are impregnated. This latter behavior - 
impregnation of the surface roughness- occurs for a few 
devices at long times and is likely due to defects at the 
edges of the channel (as suggested by direct observation). 
These few devices therefore offer the opportunity to also 
study slippage in the Wenzel state. When the CNT sur- 
face is in the Cassie state, the liquid interface lies at the 
top of the nanotube bundles, to which adsorbed beads 
are therefore restricted, fixing thereby the effective solid- 
liquid interface. On the other hand, in the Wenzel state 
beads are distributed between the top and bottom of the 
nanotubes and the "effective" (solid-liquid) interface is 
in that case fixed at the mean position between the top 
and the bottom of the CNT. The CNTs present a length 
polydispersity, responsible for the uncertainty in the lo- 
cation of the liquid interface. This uncertainty is mea- 
sured on each CNT surface and is typically of the order of 
zb200nm, similar to the characteristic deviation in height 
of the grown CNT forests (see figs El and . 

Fig. ISl shows a zoom of the water velocity profile mea- 
sured in the vicinity of the liquid interface when the CNT 
surface is in the Cassie state. The two vertical dashed 
lines represent the uncertainty in the position of the liq- 
uid interface, as discussed above. As shown in this figure, 
the velocity profile does not vanish at the boundary and 
does exhibit slippage at the interface. This slippage can 
be accounted for by the standard Navier BC for the veloc- 
ity field at the interface bdzV = v: bis the slip length and 
can be interpreted as the extrapolated distance at which 
the velocity would vanish. Experimentally, b is obtained 
by fitting the velocity profile using the Poiseuille pre- 
diction. For each surface, this procedure is repeated on 
about 5 measurement windows chosen in different loca- 
tions at the CNT surface in the microchannel from which 
an averaged value and error bar are deduced for the slip 
length. 

As mentioned above, the characteristics of the under- 
lying CNT surfaces can be varied: the distance between 
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FIG. 3: Zoom of the velocity profile in the region close 
to the CNT surfaces (in the Cassie state), as measured by 
the /i-PIV method (o), normalized by the velocity in the 
center of the channel vq (vq — 360/im/s). The errorbar 
on the velocity is ^ lO/xm.s"^, of the order of the size 
of the symbols (see 20]). The two vertical dashed lines 
correspond to the extrema of the measured position of the 
CNT surface, while the vertical solid line is the averaged 
position of the interface. The dashed-dotted line is the fitted 
velocity profile, using the parabolic (Poiseuille) prediction 
and the Navier BC. Inset: the three different profiles cor- 
respond to CNT surfaces with different roughness length 
scale L. From bottom to top, L — 1.7, 3.5, 6 /xm, and 
vo = 550; 760; 360/im/s. Slippage is shown to increase with L. 

bundles is mainly tuned by the length and density of 
the forest before functionalization in the liquid phase 
|r^ . We have examined the flow behavior on surfaces 
with various lateral length scales L, measured by com- 
puting the intensity autocorrelation of the top images of 
the SEM pictures. In the examined surfaces this typi- 
cal lateral characteristic length of the roughness varied 
between 1/im up to 6/im. We have therefore measured 
the slip length on these surfaces using the procedure de- 
scribed above. Resulting flow profiles for the Cassie state 
are shown in the inset of fig. 01 The slippage increases 
with the roughness length scale L as is summarized in 
fig. Elwhere the measured average slip length b is plotted 
against L, the error bar corresponds to the dispersion 
in the measurements. We have also measured the slip 
length for different L for systems in the Wenzel state. As 
shown in fig. El slippage is strongly reduced in the Wen- 
zel state, due to the increased liquid-solid friction in the 
absence of trapped air in the CNT, in agreement with 
theoretical predictions |T]| . 

Our results point out the absence of slippage in the 
Wenzel sate and the direct, linear, dependence of the 
measured slip length on the CNT lateral roughness scale 
L in the Cassie superhydrophobic state. This confirms 
recent theoretical predictions made for the effective slip 
length on sup erhydrophobic micro-structured surfaces 
0, El E3- In the Cassie state, the liquid is in di- 
rect contact with the solid only through a small frac- 
tion {(j)s) of the surface. At the hydrodynamic level. 



4 



1.5 



I 1 



0.5 







4' 


A 


\ + + 







2 4 6 

L (jim) 

FIG. 4: Measured slip length on the CNT surfaces as a func- 
tion of the characteristic lateral length scale L of the underly- 
ing CNT forest, o : surfaces in the Cassie state; ■ : surfaces 
in the Wenzel state; ▲ : benchmark measurement on the bare 
silicium (without CNT forest). The dotted hne is the the- 
oretical prediction for the shp length, beff = a{(l)s)L, with 
a — 0.28 (computed from the best Unear fit), 
this composite surface leads to a spatially dependent BC 
with a no-slip BC on the fraction of real solid-liquid 
interface (with fraction (j)s)^ while the remaining sur- 
face is characterized by a perfect slip BC {h = oo) on 
the suspended liquid- vapor interface. At a qualitative 
level, the theoretical results can be summarized by stat- 
ing that the effective slip length 6eff essentially saturates 
at the value fixed by the lateral scale of the roughness 
L, 6eff = Oi{(j)s)L. For instance a flow perpendicular to 
stripes yields <^ = ^ lc)g(l/ cos[^(l — ^5)]) show- 
ing a weak dependency on (ps for feasible surfaces. In 
our case, the experimental surface is closer to an array 
of pillars. Since no analytical form is available in this 
case, we have extended our previous theoretical analy- 
sis of Ref. ^3 ^o this geometry, sho wing that the main 
features of the effects described in [ifll, (for a sur- 
face made of stripes) are preserved for a pattern of pillars 
In particular, our predictions in this geometry are 
again summarized by 6eff = Oi{(j)s)L^ with a slowly vary- 
ing a{(j)s)i except very close to (ps = 0- This prediction is 
shown in fig. ^to be in good agreement with experimen- 
tal results, with a value of a :^ 0.28, corresponding to a 
constant value for (j)s of 0.15-0.2. The latter is reason- 
ably consistent with the experimental value of (j)s ^ 0.1 
estimated from the contact angle on the CNT surfaces 
165°) and the Cassie theory [l^. The slightly higher 
value of ^5 obtained in the slippage measurements may 
result from the random distribution of pillars, as opposed 
to the regular one used here in the theoretical description 
of slippage. 

Finally the slip lengths we measure are substantially 
smaller than the ones reported recently in 14] using 
global rheological measurements (6eff ^ 20/im for wa- 
ter). However, unlike the present /i-PIV method capable 



of submicrometric resolution |20|, the rheological setup 
lacks the sensitivity required to measure surface effects 
characterized by slippage in the micron range 0. Our 
results moreover points out that the effect of the under- 
lying roughness can not be overlooked in determining the 
slip length on such surfaces, therefore demonstrating that 
the gas cushion model used in Ref. fails to interprete 
slippage on nanopatterned surfaces. 

To conclude we have measured the surface hydrody- 
namic properties of superhydrophobic surfaces made of 
carbon nanotube carpets incorporated into microchan- 
nels. Such surfaces are shown to exhibit a no slip bound- 
ary condition in the Wenzel state and slip lengths of a few 
microns in the superhydrophobic (Cassie) state. More- 
over, these slip lengths were found to vary linearly with 
the lateral roughness scale L in agreement with theoret- 
ical predictions. For a microsystem with given experi- 
mental constraints (working pressure, dimensions), this 
now opens the possibility of designing optimized surfaces 
in terms of friction reduction and robustness. 

We thank K. Stephan R. Ferrigno, S. Purcell and C. 
Barentin for their help during this work. 



[1] 
[2] 

[3] 
[4] 

[s: 

[6] 
[7] 

[10: 
[11 
[12 
[13 
[14 
[15 

[16 
[17 

[is: 

[19 
[20; 



Phys. Rev. Lett 96, 
Bocquet, E. Charlaix, 
Phys. Fluids 16, 4635 



Electronic address: 'lyderic.bocquet@univ-lyonl.fr' 
T. Squires, S. Quake, Rev. Mod. Phys. 77, 977 (2005). 
E. Lauga, M. Brenner, H. Stone, Handbook of Experi- 
mental Fluid Dynamics (Springer, 2005). 
J.L. Barrat, L. Bocquet, Phys. Rev. Lett. 82, 4671 (1999) 
C. Cottin-Bizonne, B. Cross, A. Steinberger and E. Char- 
laix,P%5. Rev. Lett, 94, 056102 (2005). 
O.I. Vinogradova, G.E. Yakubov, Langmuir 19, 1227 
(2003) 

L. Joly, C. Ybert, L. Bocquet, 
046101 (2006). 

C. Cottin-Bizonne, J.-L. Barrat, L 
Nature Materials 2, 237 (2003). 
J. Ou, J.B. Perot, J.P. Rothstein, 
(2005); ibid. 17, 103606 (2005). 
E. Lauga, H. Stone, J. Fluid. Mech. 489, 55 (2003). 
J.R. Philip, Z.A.M.P. 23, 353 (1972). 
C. Cottin-Bizonne et al. Eur. Phys. J. E 15, 427 (2004). 

C. Ybert et a/., in preparation. 

D. Quere, Nature Materials 1, 14-15 (2002). 

C.-H. Choi, C.J. Kim, Phys. Rev. Lett 96, 066001 (2006). 
L. Bocquet, P. Tabeling, S. Mannevihe, to appear in 
Phys. Rev. Lett (2006). 

K.K.S. Lau et al., Nano Lett 3, (12) 1701 (2003). 

C. Journet et a/., Europhys. Lett 71, 104 (2005). 

J. Bico, B. Roman, L. Mouhn, A. Boudaoud, Nature 432, 

690 (2004). 

P. Vulto et al., Lab on a Chip 5, 158 (2005). 

P. Joseph, P. Tabeling, Phys. Rev. E, 71, 035303 (2005). 



